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Abstract 

We perform a detailed study of the Earth matter effects on supernova neutrinos 
with neutrino oscillation parameter LMA and small #13. The Earth effects show 
significant dependences on the neutrino path length inside the Earth and the value 
of 

A77t^2- We investigate rather optimistically a possibility that we can probe the 
value of Awij 2 by the Earth effects. We assume that 812 and the direction of the 
supernova are known with enough accuracy and that the resonance that occurs 
at higher density in supernova envelope is completely nonadiabatic. Further the 
neutrino spectra before neutrinos go through the Earth are assumed to be known. 
Then we show that making use of these dependences, we can obtain implication 
for the value of Amf 2 by comparing the observed energy spectrum to the predicted 
one. When SK detects neutrinos from supernova at lOkpc which traveled through 
the Earth (nadir angle < 80 degree), Amf 2 can be determined with an accuracy 
of ~ 10%. In much of the neutrino-detection-time-Amf 2 plane, AmJ 2 might be 
determined with an accuracy equal to or better than ±0.5 x 10 _5 eV 2 . 
PACS : 14.60.Pq; 14.60.Lm; 96.40.Tv; 97.60. Bw; 
Keywords : Neutrino oscillation; Supernovae; Earth effects; 
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1 Introduction 



A type II supernova is a prodigious source of neutrinos. Almost all of the binding energy, 
Ef, ~ 10 53 erg, is radiated away as neutrinos. These neutrinos carry information about 
both the core collapse process and intrinsic properties of the neutrinos. In fact, the 
detection of the neutrino burst from SN1987A |1|, Q induced a number of studies on these 
areas |, |,|,§. 

Neutrino oscillation, which is confirmed by the observations of solar and atmospheric 
neutrinos, can affect the energy spectrum of supernova neutrino drastically. Neutrinos 
of all flavors are produced in the high dense region of the iron core |7J and interact with 
matter before emerging from the supernova. The presence of non-zero masses and mixing 
in vacuum among various neutrino flavors results in strong matter dependent effects, 
including conversion from one flavor to another. Hence, the observed neutrino flux in 
the detectors may be dramatically different for certain neutrino flavors and for certain 
values of mixing parameters, due to neutrino oscillation. Similar matter effect occurs in 
the Earth, too. 

These oscillation effects on neutrino spectra depend strongly on neutrino oscillation 
parameters: mixing angles and mass spectrum. Some of them are firmly established but 
others are not. All the existing experimental results on the atmospheric neutrinos can be 
well described in terms of the <-> v T vacuum oscillation with mass squared difference 
and the mixing angle given by ||, 

A^atm ~ 7 x l(r 3 eV 2 , sin 2 29 « 1. (1) 

In contrast, for the observed v e suppression of solar neutrinos four solutions are still 
possible [[|, 110, O]: large mixing angle (LMA), small mixing angle (SMA), low Am 2 



(LOW), and vacuum oscillation (VO). 



(LMA) Am| « (1 ~ 10) x 10" 5 eV 2 , sin 2 20 & « 0.7 ~ 0.95 (2) 

(SMA) Am% » (4 ~ 10) x 10~ 6 eV 2 , sin 2 26 & « (2 ~ 10) x 10~ 3 (3) 

(LOW) Am 2 Q « (0.5 ~ 2) x 10~ 7 eV 2 , sin 2 2# « 0.9 ~ 1.0 (4) 
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(VO) Ami « (0.6 ~ 6) x 10~ 10 eV 2 , sin 2 26 Q « 0.8 ~ 1.0 (5) 



For #i3, the mixing angle between mass eigenstate Ui, 1/3, only upper bound has been known 
from reactor experiment and combined three generation analysis |13|, [14], [15| . Also the 



nature of neutrino mass hierarchy (normal or inverted) is still a matter of controversy [16 

Supernova is a completely different system from solar, atmospheric, accelerator, and 
reactor neutrinos in regard to neutrino energy and flavor of produced neutrinos, propa- 
gation length and so forth. Then neutrino emission from a supernova is expected to give 
valuable information that can not be obtained from neutrinos from other sources. 

There have been some studies on the supernova neutrino as an oscillation parameter 
prober. Dighe and Smirnov |TBJ have studied the role the supernova neutrinos can play in 



the reconstruction of the neutrino mass spectrum. Dutta et al. [19| showed numerically 
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that the events involving oxygen targets increase dramatically when there is neutrino 
mixing. We have shown in our previous work p(| that the degeneracy of the solutions 
of the solar neutrino problem can be broken by the combination of the SK and SNO 
detections of a future Galactic supernova. 



The Earth effects have also been studied by several authors. In our previous work [21 



we analyzed numerically the time-integrated energy spectra in a simple case and found the 
possibility to probe the mixing angle #13. Lunardini and Smirnov p3] performed a study 
of the Earth effects taking the arrival directions of neutrinos into account and showed that 
studies of the Earth effects will select the solution of the solar neutrino problem, probe 
the mixing U e s and identify the hierarchy of the neutrino mass spectrum. 

In this paper, we perform a detailed study of the Earth matter effects on supernova 
neutrinos and show that the detection of Earth matter effects allows us to probe Am^ 
more accurately than by solar neutrino observations. 



2 Three-flavor formulation 

In the framework of three-flavor neutrino oscillation, the time evolution equation of the 
neutrino wave functions can be written as follows: 
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where A(t) = y/2Gpn e (t), Gp is Fermi constant, n e (t) is the electron number density, 
Am| is the mass squared differences, and E is the neutrino energy. In case of antineutrino, 
the sign of A(t) changes. Here U is a unitary 3x3 mixing matrix in vacuum: 



U 



C12C13 S12C13 S13 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 
S12S23 - C12C23S13 -C12S23 _ S12C23S13 C23C13 



(8) 



where s» 



sin 6ij , Cij 



cos Qi 



l,2,3(i < j). We have here put the CP phase 



, t j — o±±.L , — ^vjowjj for i, j 
equal to zero in the CKM matrix. 

In H(t), the first term is the origin of vacuum oscillation, and the second term A(t), 
which is the only time-dependent term in H(t), is the origin of MSW effect. 



3 Determination of neutrino oscillation parameters 

In this section we summarize our previous studies on the effects of neutrino oscillation on 
supernova neutrino and show that future detection of neutrino from Galactic supernova 
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allows us to break the degeneracy of the solutions of the solar neutrino problem and to 
probe #13. 

In three- flavor neutrino oscillation in the star is studied. We calculated the ex- 
pected event rate and energy spectra, and their time evolution at the SuperKamiokande(SK) 
and the Sudbury Neutrino Observatory(SNO). For the calculations of them, we used a 
realistic neutrino burst model based on numerical simulations of supernova explosions 
|7|, |23j and employed a realistic density profile based on a presupernova model |2|| . The 



distance between the supernova and the Earth was set to lOkpc. 

By solving numerically the differential equations (^) from the center of supernova to 
the outside of supernova, we obtained conversion probabilities P{v a ^p), i.e., probabilities 
that a neutrino of flavor a produced at the center of supernova is observed as a neutrino 
of flavor (3. 

We assumed normal mass hierarchy and used four sets of mixing parameters shown in 
table |[ Here # 12 and Am^ correspond to the solutions of solar neutrino problem and #23 
and Am^ 3 correspond to the solution of atmospheric neutrino problem. The value of #13 



is taken to be consistent with current upper bound from reactor experiment ||12|| . These 
models are named after their values of mixing angle: LMA-L means that #12 is set to be 
LMA of solar neutrino problem and #13 is large. 

We then found that when there is neutrino oscillation, neutrino spectra are harder 
than those in absence of neutrino oscillation. This is because average energies of v e and 
v e are smaller than those of v x (either of u^, z/ r , and their antineutrinos), and neutrino 
oscillation produces high energy v e and v e which was originally v x . This feature was used 
as a criterion of magnitude of neutrino oscillation. We calculated the following ratio of 
events at both detectors: 

„ number of events at 30 < E < 70MeV , , 

Rsk = (9) 

number of events at 5 < E < 20MeV 

^ _ number of events at 25 < E < 70MeV 
number of events at 5 < E < 20MeV 
The plots of Rsk vs. -Rsno are shown in Fig.|l| The errorbars include only statistical 
errors. The difference among the following three groups is clear: (l)LMA-L and LMA-S, 
(2)SMA-L, and (3)SMA-S and no oscillation. 



In |21|, the Earth effects on the energy spectra are studied using the result of |20 
For mass squared differences in table |], the supernova neutrinos arrive at Earth in mass 
eigenstates. This is because neutrino eigenstates with different masses lose coherence on 
the way from the supernova to the Earth. Since the eigenstates of Hamiltonian in matter 
differ from the mass eigenstates in vacuum, supernova neutrinos begin to oscillate again 
in the Earth. The numbers of neutrinos of each mass eigenstates at the surface of the 



Earth are determined by neutrino oscillation in the supernova |2T . 

We analyzed numerically the time-integrated energy spectra of neutrino in a mantle- 
core-mantle step function model of the Earth's matter density profile. We assumed that 
neutrino arrived at the detectors after traveling through the Earth along its diameter. We 
then found that the Earth matter affects the v e spectrum significantly in model LMA-S 
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and u e spectrum slightly in model LMA-S and LMA-L. In other models, there are no 
significant Earth effects. We concluded that we can differentiate LMA-L from LMA-S, by 
observing the Earth matter effect in v e . 

Thus we can distinguish all the four models in table |]. But there are some ambiguities 
besides statistical errors. One is the mass of the progenitor star. Supernovae with different 
progenitor masses may result in different original neutrino spectra and neutrino oscillation 
effects. Studies on this point are now in progress. But dependence of shape of neutrino 



spectra on progenitor mass is not so large |25| and we would still be able to distinguish 
the models. 

Supernova model which we use does not consider rotation of the progenitor star. In 
general, however, stars rotate through their lives and recent numerical simulation indicates 
that the rotation facilitate supernova explosion PB| . Rotation of the progenitor star can 
affect its density profile and the dynamics of the neutrino oscillation in supernova. 

Another ambiguity is the direction of supernova. The trajectories of neutrinos change 
according to the location of the supernova, the position of the detector and the time t of 
the day at which the burst arrives at the Earth [^| . Fig shows the dependence of the 
length d of neutrino trajectory on the time t for the three detectors, SK, SNO and LVD. 
d is in unit of the diameter of the Earth. Neutrino trajectory can also be described by 
the nadir angle 9 n of the supernova with respect to the detector. In the next section, we 
perform a detailed study of the Earth matter effect. 

4 Detailed study of the Earth matter effect 

In this section, we concentrate on model LMA-S and discuss the dependence of the Earth 
matter effect on the distance which neutrinos travel in the Earth and Am^. We consider 
three detectors: SK, SNO and LVD. But energy spectra at SK and LVD in model LMA-L 
are the same as in model LMA-S, since most of the events at SK and LVD are induced by 
z/ e s which are influenced by the Earth matter in both LMA-L and LMA-S. The dependence 
on Amf 2 is expected because the oscillation length in the Earth between mass eigenstates 



i>i and V2 is comparable (~ 1000km) to the diameter of the Earth in LMA-S model [|21 



Here we use the realistic Earth density profile |27j , while we used step-like model in our 



previous work |[21|| . By solving (|j) along the Earth density profile numerically, we obtain 



conversion probability P Vi ^ Va , i.e. probability that ith mass eigenstate at the surface of 
the Earth is detected as neutrino of flavor a. Then combining these probabilities with 



the result of neutrino oscillation in supernova f20 |, we obtain neutrino flux and at each 
detector. 

4.1 neutrino events at SuperKamiokande 

SuperKamiokande (SK) is a water Cherenkov detector with 32,000 ton pure water based 
at Kamioka in Japan. The relevant interactions of neutrinos with water are as follows: 

v e +p -> n + e + (CC) (11) 
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(12) 
(13) 
(14) 
(15) 
(16) 

where CC and NC stand for charged current and neutral current interactions, respectively. 
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For the cross sections of these interactions, we refer to |28[ . 

SK is not available now due to the unfortunate accident |29| but is expected to be 
restored within the year. Since the number of PMTs will decrease to half the original, 
energy threshold and energy resolution will become worse. The detection efficiency is 



expected to be 90% at 8MeV|30|]. In the above interactions, the v e p CC interaction 
[Eq.flTTDl has the largest contribution to the detected events at SK. Hence the energy 
spectrum detected at SK (including all the reactions) is almost the same as the spectrum 
derived from the interaction Eq.flTTD only. 

4.2 neutrino events at SNO 

Sudbury Neutrino Observatory(SNO) is a water Cherenkov detector based at Sudbury, 
Ontario. SNO is unique in its use of 1000 tons of heavy water, by which both the charged- 
current and neutral-current interactions can be detected. The interactions of neutrinos 
with heavy water are as follows, 

v e + d -> p + p + e~ (CC) (17) 

u e + d -> n + n + e + (CC) (18) 

v x + d -> n + p + u x (NC) (19) 

i7 x + d -> n + p + v x (NC) (20) 

The two interactions written in Eqs.(pj]) and ( [Tq) are detected when electrons emit 
Cherenkov light. These reactions produce electrons and positrons whose energies sen- 
sitive to the neutrino energy, and hence the energy spectra of electrons and positrons give 
us the information on the original neutrino flux. In this work, we mainly take into account 
these two charged current interactions, since the number of neutral current events does 



not change by neutrino oscillation. For the cross sections, we refer to [^T||. The efficiency 
of detection is set to be one, because we have no information about it. 

The SNO detector has also 7,000 tons of light water which can be used to detect 
neutrinos. This can be considered to be a miniature of SuperKamiokande (32,000 tons of 
light water). Then the number of events detected by light water at SNO is 7/32 of that 
at SuperKamiokande. 

4.3 neutrino events at LVD 

The Large Volume Detector (LVD) in the Gran Sasso Underground Laboratory is a v 
observatory mainly designed to study low energy neutrinos from gravitational stellar col- 



6 



lapse. LVD consists of 2 kinds of detectors, namely: liquid scintillator, for a total mass 
of 1840 tons, and streamer tubes for a total surface of about 7000m 2 . 
The bulk of events is due to the capture reaction: 

v e +p-*n + e + n + p -> d + 7 2 . 2 MeV (21) 

Further, about 5% of the events are due to neutral current interactions with 12 C which 
deexcites emitting a 15.1 MeV 7. Moreover, 3% of events are due to elastic scattering of 
all neutrino flavours on electrons, and less than 1% to charged current interactions of u e 
and v e with 12 C nuclei. For the cross sections, we refer to [ 32| . The appropriate detection 



efficiency curve is also taken into account [22 . 



4.4 nadir angle dependence of energy spectra 

Fig.^|, [| and || shows the nadir angle dependence of neutrino spectrum at SK, SNO and 
LVD, respectively. Neutrino oscillation parameters are set to the values in model LMA-S 
except Am 2 2 = 2 x 10 _5 eV 2 . Since the oscillation lengths are comparable to the radius 
of the Earth, the Earth effects are strongly dependent on the nadir angle. 



5 AmJ 2 dependence of energy spectra and determi- 
nation of Amf 2 

5.1 method 

Fig.|| [7| and [8] shows the Am 2 2 dependence of energy spectrum at SK, SNO and LVD. 
Neutrino oscillation parameters are set to the values in model LMA-S except Am 2 2 . Nadir 
angle is set to degree. The larger Am 2 2 results in higher-frequency oscillation in the 
energy spectra with respect to the energy. This is because the neutrino oscillation length 
is proportional to the inverse of Am 2 . 

As can be seen, there is significant dependence of the Earth effect on Am 2 2 . Making 
use of this dependence, we can probe Am 2 2 by the observations of the supernova neutrino 
spectra more accurately than by the observations of solar neutrino, if neutrino oscillation 
parameters are as in model LMA-S. 

We assume that the direction of the supernova is given by direct optical observations 
or by the experimental study of the neutrino scattering on electrons |33|, [34|]. Other 
parameters of neutrino oscillation except Am 2 2 are also assumed to be known with enough 
accuracy. But in fact, #23 does not affect the Earth matter effects and #13 affects only 
the magnitude of the Earth matter effects through the adiabaticity of the H-resonance in 
supernova envelope |8j. If sin 2 13 < 10" 5 , the H -resonance is completely nonadiabatic for 



broad range of Am| 3 and in this case the Earth matter effects do not depend on the value 
of #13. So what we need with enough accuracy is 9 12 only if we assume the H-resonance 
is completely nonadiabatic. 

Comparing the observed neutrino spectra to the predicted spectra of various values 
of Am 2 2 , we can determine the value of Am 2 2 . We perform Monte Carlo simulation to 
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obtain the expected observed spectra at each detector. Here we take statistical fluctuation 
in number of neutrino events and energy resolution of each detector into account. Energy 
spectra are then binned according to the energy resolution. For the value of the energy 
resolution, we refer to p2| but that of SK after restoration is expected to be a/2 times as 



bad as before the accident pCf . The energy resolution of each detector is shown in Fig.g. 

The simulated spectrum is compared to the predicted spectrum with a x 2 method. 
The definition of the reduced \ 2 ■, which is a function with respect to Am 



2 

12 > 



Y 2 fAmM--V^ 1 1 12,) (22) 

where d is the number of bins and Nf ira and Nf re (Am\ 2 ) are the simulated and predicted 
number of events in zth bin, respectively. The Am\ 2 with the least x 2 {Am\ 2 ) is considered 
to be the value determined by this method. 

5.2 representative results 

We make 1000 "supernovae" and "observe" the neutrino spectra assuming that the true 
value of Amf 2 is 4 x 10 _5 eV 2 . Fig.[T(] show the representative results of Am\ 2 determi- 



nation by the method described above. These are relative frequencies that the observed 
spectrum is identified with the theoretical spectrum with the various values of Am\ 2 . The 
nadir angle is set to be degree at each detector and the \ 2 method is performed with 
the data of each detector only. 

As can be seen, if the supernova occur at the direction which the nadir angle is degree 
at SK or SNO, Am\ 2 can be determined with a high accuracy (4.00 ± 0.07) x 10~ 5 eV 2 or 
(4.00 ± 0.20) x 10~ 5 eV 2 in case of SK and SNO, respectively. However, in case of LVD, it 
is difficult to determine Am\ 2 . The difference in the accuracy of determination of Am\ 2 
is due to the difference in the properties of the detectors. Here the properties of the 
detectors mean the energy resolution, the detectability of v e and the number of events. 

Since the Earth effects are large in v e channel and relatively small in v e channel, SK 
and LVD, at which most of events are i> e , are disadvantageous. But the statistical error 
is very small at SK due to the large event number, ~ 10000 (10 kpc), and SK is the most 
efficient to identify the Earth effects. SNO is also efficient due to its high energy-resolution 
and detectability of u e although the number of events is small. Then LVD is not favorable 
for detecting the Earth effects. The properties of the detectors are summarized in the 
Table |. 

It should be noted that the effectiveness of determining Am\ 2 depend on the value of 
Am\ 2 itself and nadir angle. These dependences are discussed in the following subsection. 



5.3 determination of Am? 
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Figjfl] shows the relative frequencies at SNO that the simulated spectrum based on the 
theoretical spectrum with Am\ 2 = 2(solid), 4(dashed), 6(long — dashed) x 10~ 5 eV 2 is iden- 
tified with the theoretical spectrum with the various values of Am\ 2 . The nadir angle is 
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set to be degree. When Am 2 2 is large, the oscillation in energy spectrum becomes rapid 
and then statistical errors make it difficult to identify the Earth effect and determine 
Am\ 2 . However at SK, its large number of events allows us to determine Am 2 2 rather 
accurately irrespective of the true value of Am 2 2 (Fig.|i~2|). 

Fig.[13| shows the relative frequencies at SNO at various nadir angles: 15 (solid), 80 
(dashed) and 85 (long-dashed). The value of Am 2 2 is set to be 4 x 10~ 5 eV 2 . Since the 
oscillation length of 60 MeV v in Earth matter is ~ 6000km |^T[, the Earth effects do not 
appear unless nadir angle is smaller than ~ 80 degree. 

Now we consider the problem, 'How accurately can we determine Amj 2 when the true 
value of AmJ 2 and the arrival time of neutrinos are given?' 

As is discussed above, when the nadir angle at SK is smaller than 80 degree (time 
= ~ 8 and 19 ~ 24 hour), the value of Am^ can be determined accurately using only 
the SK data irrespective of the value of Am 2 2 in the allowed region of the solar neutrino 
problem. When the nadir angle is smaller than 80 degree at SNO and greater than 80 
degree at SK (time = 10 ~ 19 hour), the value of Am 2 2 can be determined accurately 
using only the SNO data if the true value of Am 2 2 is smaller than ~ 4x 10~ 5 eV 2 . However, 
if Am 2 12 > 4 x 10~ 5 eV 2 , the single data of SNO is not sufficient. We performed then a x 
method combining data of various detectors for time = 8 ~ 19 hour. 

Fig.[14] is the contour map of the probability that the value Am 2 2 can be determined 
with an accuracy equal to or better than ±0.5 x 10 _5 eV 2 . The probability is more than 
0.7 for three fifths of the time-Am 2 2 plane. 



6 Discussion and conclusion 

In this paper, we performed a detailed study of the Earth matter effects on supernova 
neutrinos with neutrino oscillation parameter LMA and small # 13 . We showed that we 
can probe Am 2 2 accurately by comparing the observed energy spectra to the predicted 
one. In much of the time-Am 2 2 plane, Am 2 2 can be determined with an accuracy equal 
to or better than ±0.5 x 10~ 5 eV 2 . When SK detect neutrinos from supernova at lOkpc 
which traveled through the Earth (nadir angle < 80 degree), Am 2 2 can be determined 
with an accuracy of ~ 10%. This accuracy is amazing compared to that from current 
solar neutrino experiments, which determine only the order of Am 2 2 . But as we stated 
before, # 12 and the direction of the supernova need to be known with enough accuracy. 
How much our results are affected by the uncertainties of the other parameters will be 
our future work. 

Am 2 2 can also be probed by KamLAND experiment [|35fl . If the solution of solar 
neutrino problem is LMA, Am 2 2 is expected to be known by three-year data-taking as 
precisely as by supernova neutrino discussed above. There is some possibility that we can 
know it from a Galactic supernova before the completion of the KamLAND experiment. 
Even if the results of KamLAND experiment is earlier than future galactic supernova, 
it is important to cross-check the results of ground-based experiments by astrophysical 
observations. Furthermore, it is meaningful itself to notice that we can obtain implication 
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for Am 2 from supernova neutrino. 

We assumed in this paper that the neutrino spectra at the surface of the Earth were 
given. But since the original fluxes based on the numerical supernova model have some 
ambiguities, the sensitivity to Am] 2 might be worse than our estimation. Estimation 
of uncertainties of the numerical supernova model is hard because there is currently no 



successful simulation other than that by the Lawrence Livermore group [23]. There are, 
however, some studies about the temperatures of the produced neutrinos. According to 
them, they are typically [|36{|. 

(E Ue ) = 10 - 12MeV, (E Pe ) = 14 - 17MeV, (E Vx ) = 24 - 27MeV. (23) 

On the other hand, in the model that we used, they are 

(E Ve ) ~ 13MeV, (E Pe ) ~ 16MeV, (E Vx ) ~ 23MeV, (24) 

although the spectra are not exactly the Fermi-Dirac distribution ]7|. These values seem 
close to the typical values and the analysis taking these uncertainties into account will be 
performed in our future work. 

Furthermore, supernova neutrino spectra depend on the mass of the progenitor star. 
However, we will be able to reconstruct the neutrino fluxes at the surface of the Earth 
from the data from the detector which detect the neutrinos directly from the supernova. 
Although the statistical errors then become larger and the accuracy of the determination 
of Amf 2 becomes worse, the Earth effects on the supernova neutrinos will still give valuable 
information which cannot be obtained from the other neutrino sources. 
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Figure 1: The plot of Rsk vs. Rsno for all the models 
the statistical errors. 
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Figure 2: The distance which neutrinos travel in the Earth as functions of the arrival 



time of neutrino burst |22| . We assume a supernova which located in the Galactic center. 
We fixed t = as the time at which the star is aligned with the Greenwich meridian. 
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Figure 3: Nadir angle dependence of energy spectrum at SuperKamiokande. Neutrino 
oscillation parameters are set to the values in model LMA-S except Amf 2 = 2 x 10~ 5 eV 2 . 
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Figure 4: Nadir angle dependence of energy spectrum at SNO taking only CC events into 
account. Neutrino oscillation parameters are set to the values in model LMA-S except 
Am? 9 = 2 x 10" 5 eV 2 . 
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Figure 5: Nadir angle dependence of energy spectrum at LVD. Neutrino oscillation 
parameters are set to the values in model LMA-S except Amf 2 = 2 x 10~ 5 eV 2 . 
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Figure 6: Am^ dependence of energy spectrum at SuperKamiokande. Neutrino oscilla- 
tion parameters are set to the values in model LMA-S except Am^. Nadir angle is set to 
degree. 
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Figure 7: Amj 2 dependence of energy spectrum at SNO taking only CC events into 
account. Neutrino oscillation parameters are set to the values in model LMA-S except 
Aml 2 . Nadir angle is set to degree. 
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Figure 8: Am^ dependence of energy spectrum at LVD. Neutrino oscillation parameters 
are set to the values in model LMA-S except Am^. Nadir angle is set to degree. 
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Figure 9: Energy resolution at each detector: SK (solid), SNO (dotted) and LVD 
(dashed) |22[. 
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Figure 10: Relative frequencies that the simulated spectrum based on the theoretical 
spectrum with Am^ = 4 x 10 _5 eV 2 is identified with the theoretical spectrum with the 
various values of Am^ 2 . The nadir angle is set to be degree at each detector and the x 2 
method is performed with the data of each detector only. 
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Figure 11: Relative frequencies at SNO that the simulated spectrum based on the theo- 
retical spectrum with AmJ 2 = 2 (solid), 4 (dashed), 6 (long — dashed) x lCT 5 eV 2 is identified 
with the theoretical spectrum with the various values of Am^. The nadir angle is set to 
be 80 degree. 
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Figure 12: The same as in Fig. 11 for SK 
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Figure 13: Relative frequencies at SNO at various nadir angles: 15 (solid), 80 (dashed) 
and 85 (long-dashed). The value of Am 2 2 is set to be 4 x 10~ 5 eV 2 . 
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Figure 14: Contour map of the probability that the value Amf 2 can be determined with 
an accuracy equal to or better than ±0.5 x 10~ 5 eV 2 . 
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Table 1: Sets of mixing parameter for calculation 



model 


sin 2 29 12 


sin 2 29 23 


sin 2 29 13 


LMA-L 


0.87 


1.0 


0.043 


LMA-S 


0.87 


1.0 


1.0 x 10" 


SMA-L 


5.0 x 10" 3 


1.0 


0.043 


SMA-S 


5.0 x 10~ 3 


1.0 


1.0 x 10" 



7.0 x 10" 5 3.2 x 10" 3 ' LMA 

7.0 x 10~ 5 3.2 x 10" 3 LMA 

6.0 x 10~ 6 3.2 x 10" 3 SMA 

6.0 x 10~ 6 3.2 x 10" 3 SMA 



Table 2: Properties of the detectors 

detector SNO SK LVD 

main event u e , v e v e v e 

number of events (lOkpc) 300 10000 800 

energy resolution < 3 MeV < 6 MeV < 7 MeV 
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